Shape Dynamics of Interfacial Front in Rotating Cylinders 
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The evolution of the interface propagation in a slowly ro- 
tating half-filled horizontal cylinder is studied using MRI. Ini- 
tially, the cylinder contains two axially segregated bands of 
small and large particles with a sharp interface. The pro- 
cess of the formation of the radial core is clearly captured, 
and the shape and the velocity of the propagating front are 
calculated by assuming a one-dimensional diffusion process 
along the rotation axis of the cylinder and a separation of 
time scales associated with segregation in the radial and axial 
directions. We found that the interfacial dynamics are best 
described when a concentration dependent diffusion process 
is assumed. 
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I. INTRODUCTION 

When granular materials are put in a horizontal rotat- 
ing cylinder, most particles undergo a solid body rotation 
except for a thin layer of particles flowing near the free 
surface. With increasing rotation rate, particles first ex- 
hibit intermittent avalanches, then continuous flow near 
the surface and finally a ring structure due to centrifugal 
force In industrial applications, rotating cylinders are 
primarily used for mixing different components But 
it is well known that particles of different sizes or density 
show radial segregation on small time scales and axial 
segregation on large time scales which compete against 
the mixing process 

One of the ways to investigate the three-dimensional 
particle motion is to stop the rotating cylinder and take 
samples from different locations JlO[. However, this in- 
evitably involves a partial destruction of the internal 
structure which might provide vital information. An- 
other shortfall of this approach is that it is time consum- 
ing to obtain a global information about the particle mo- 
tion. Recently, non-invasive magnetic resonance imaging 
(MRI) has been used to interrogate the internal struc- 
tures of evolving granular assemblies |ll[]. Since the cur- 
rent MRI fiow measurement is designed for a steady flow 
of particles and recording a single image can take several 
minutes, the cylinder is also stopped to capture a series 
of static images of the propagating front starting with 
the same set of initial conditions. Nakagawa et al. 
first applied MRI to study the dynamics of granular flow 
in a rotating cylinder and recorded velocity and density 
profiles along and perpendicular to the flowing layer of 
particles. Using a spin-tagging sequence, the flow and 



diffusion in a vertically shaken container of poppy seeds 
was recently measured and convection rolls found [ p^ . 
Most recently, the shape of the radially and axially seg- 
regated core were also investigated by MRI [p|jT^]. 

Not much is known yet about how this radially seg- 
regated core evolves. In the common segregation exper- 
iments, one starts with a well mixed initial state con- 
taining small and large particles. The formation of a ra- 
dial core of small particles which extends over the whole 
length of the cylinder is already visible after a few rota- 
tions 1^. This radial core may become unstable forming 
band patterns which is called axial segregation. This pro- 
cess was best demonstrated by Hill et al. Q using MRI. 
The final number of bands, their positions and widths 
varied from experiment to experiment but Nakagawa |^ 
found the three band configuration to be stable after an 
extended number of rotations. These bands may not be 
pure and a radial core might still be present whereas 
Chicharro et al. |Q rotated two sizes of Ottawa sand for 
two weeks at 45 rpm and found a final state of two pure 
bands each filling approximately half of the cylinder, i.e. 
no radial core was found. Recently, experimental evi- 
dences were given for two new mechanisms of axial seg- 
regation, namely traveling waves ^ and avalanche medi- 
ated transport |14| which were only observed when non- 
spherical particles were used. 

Since the dynamics of the axial segregation process de- 
pends on the homogeneity of the initial mixture, it is 
desirable to start with a better defined initial packing. 
We propose to start with a fully axially segregated state 
of a binary particle mixture and use MRI to investigate 
the three-dimensional segregation front which is initially 
seen as a sharp flat interface. We study the propagation 
velocities and compare them with results obtained from a 
simple diffusive process where a concentration dependent 
diffusion coefficient is used. The shape of the segregation 
front can be calculated approximately analytically when 
a separation of time scales for the segregation in the ra- 
dial direction and the diffusion in the axial direction is 
assumed. 




FIG. 1. A simplified sketch of the initial configuration: 
large particles are all in the left half of the cylinder and shown 
in gray. This corresponds to the first side view image in Fig. 



II. EXPERIMENTAL SETUP AND STUDIES 

The initial state was prepared in such a way that the 
left half of the cylinder contained only large particles hav- 
ing a diameter of 4 mm and the right half contained only 
small particles having a diameter of 1 mm. It is similar 
to the one used in |l5| and is sketched in Fig. Q where 
the larger particles are shown in gray. A 50-50 mixture 
by volume was used and the total volume was chosen to 
give a roughly half-filled cylinder. 

Axial migration of the core of small particles is a three 
dimensional process and its dynamics should be moni- 
tored by a non-invasive method. The very strength of 
MRI among other existing non-invasive techniques lies 
in its ability of measuring dynamic flow properties such 
as velocities and velocity fluctuations. Measurement of 
true concentration still remains to be a challenging prob- 
lem because in granular flows concentration and particle 
agitation, characterized by the granular temperature, are 
strongly coupled which is predicted by kinetic theories. 
It is not trivial to separate causes for the lesser concen- 
tration signals. Less challenging but equally important is 
MRI's ability of probing a dense static granular system. 
In fact, the denser the system of interest is, the higher 
the quality of images become. Since the MRI technique 
used is best suited for flows in a steady state, the cylinder 
was stopped when the images were taken. 

Our MRI imager/spectrometer (Nalorac Cryogenics 
Corp) has a 1.9 Tesla superconducting magnet (Oxford) 
with a bore diameter of 31 cm. The useful space inside 
the bore is a sphere of diameter of about 8 cm after the in- 
sertion of the gradient coils and the rf probe. The acrylic 
cylinder with 10 cm and 7 cm in its length and diameter, 
respectively, was rotated at a constant rate of 11.4 rpm 
by a dc servo motor (12FG) and a controller (VXA-48- 
8-8) made by PMI inside the imaging apparatus. The 
material properties, the drum geometry and its rotation 
speed were chosen such that the drum was operating in 
the beginning of the continuous flow regime. Images were 
taken after 0, 15, 30, 45, 60, 75, 90 and 600 seconds of 
rotation. The particles used were spherical pharmaceu- 
tical pills of 1 and 4 mm diameter which contain a liquid 
core of vitamin oil. Louge et al. conducted a detailed 
binary impact experiment using these particles and esti- 
mated the normal coefficient of restitution to be about 
0.9. In the perpendicular direction, however, they have 
found that these liquid-filled particles exhibit rolling con- 
tacts with negligible compliance. 

The cylinder was put into the magnet for record- 
ing three-dimensional intensity signals (64x64x64 points) 
which took around 40 minutes each. Due to the lim- 
itations in the spatial resolution, the smaller particles 
cannot be resolved individually and appear light gray in 
the images. Since the larger particles contain more liq- 
uid they give a higher signal and show up in black or 
dark gray in the images. Despite the fact that our im- 
ages are three-dimensional, we have decided to illustrate 



the dynamics of the conical shape of the migrating front 
by showing two different cross sectional views in Fig. H, 
namely the side views on the left and the top views on 
the right. Initially (0 sec), the large and small particles 
were packed in the left and right half of the cylinder, re- 
spectively, and the heaping of the large particles at the 
left end in the top left image is due to an actual uneven 
initial packing. This quickly disappeared and was flat- 
tened to the same level as in other parts of the cylinder 
once the cylinder was in motion. 



side view 



top view 








45 sec 
















FIG. 2. Cross sectional views of the interface propagation 
for a mixture of 1 and 4 mm vitamin pills in a half-filled 
rotating cylinder. The small particles cannot be resolved in- 
dividually and the region they occupy is shown in light gray. 

The side views were obtained by superimposing the 16 
vertical planes in the central region of the cylinder and 
taking the highest intensity value of all planes. This av- 
eraging increased the contrast and we verified that the 
shape of the front did not change significantly within 
these planes. For the top view, 12 planes in the lower 



half of the cyHnder were superimposed. By this proce- 
dure, we can clearly identify the regions containing large 
particles without loosing spatial information. The time 
evolution in Fig. || shows that the cascading layer is al- 
most entirely composed of large particles whereas the 
small particles can sink into the surface layer more eas- 
ily and consequently their concentration in the layer will 
decrease. The distance of the tip of the migrating front 
of small particles and the free surface, denoted by Y in 
Fig. ^, remains more or less constant. This indicates 
that the small particles in the cascading layer are re- 
sponsible for the front advancement. Immediately after 
the cylinder starts to rotate, there is a mixture of large 
and small particles flowing in the cascading layer devel- 
oped at the interface. Due to percolation mechanism, the 
smaller particles travel under the free surface and move 
in the axial direction driven by the concentration gradi- 
ent. The smaller particles advance more easily into the 
region occupied by larger particles since there are more 
voids in the cascading region for the small particles to 
move into. We also observed that the larger particles 
traveling on the surface always reached the other end of 
the cylinder before the smaller particles did. Eventually, 
this leads to an extended radially segregated core, shown 
in Fig. ^ and also observed in other experiments |^ . 



whereas the boundary conditions read 




600 sec 




FIG. 3. Side view of the segregation front for later times 
at 90 s (left) and at 600 s (right) after start of rotation. 



III. RADIAL CORE DEVELOPMENT 

Assuming random particle motion in the axial direc- 
tion (z axis), one component systems could be well de- 
scribed by a diffusion process according to Fick's Second 
Law [|l^,|l^. The interface of a two component system 
can also be studied in this fashion and the diffusion equa- 
tion reads 



dC{z,t) _ d_ r dC{z,t) 
dt dz V dz 



(1) 



where C{z,t) and D denote the relative concentration 
by volume of the smaller particles and the corresponding 
diffusion coefficient, respectively. The initial condition 
for a cylinder with length L are 
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which states that there is zero axial flux at the boundaries 
due to the end caps. 

For a constant diffusion coefficient, Eq. (|^) can be 
solved analytically for the specified initial and boundary 
conditions and the solution reads 
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Hogg et al. |lO| solved a similar diffusion problem for a 
single component system by replacing the time derivative 
with the derivative with respect to the number of rota- 
tions thus assuming that the dynamic effects on the dif- 
fusion process was not significant within the range of the 
rotation rates studied. Nakagawaet al. investigated a 
mixing process of 50-50 mixtures of liquid filled pharma- 
ceutical pills of different sizes. After different numbers of 
revolutions they inserted 12 dividers along the rotational 
axis and measured the fraction of small particles in each 
of the 13 segments. They obtained one-dimensional con- 
centration profiles similar to the theoretical prediction 
given by Eq. (^). However, for longer rotation times, 
their profiles showed a slight asymmetry with respect to 
the position of the initial front. It should be kept in 
mind, however, that the final configuration in both our 
MRI study and Ref. |l5j was a radial segregated core of 
small particles and not a mixture as in Refs. ||l0| , p^ ." 

It was found experimentally |^,^ and numerically p^ ] 
that a nearly complete radial segregation can be achieved 
after only a few rotations. In the above, successive MRI 
pictures were separated by 15 seconds which corresponds 
to roughly 3 full rotations and each particle had a chance 
to participate in approximately 6 avalanches. In order to 
calculate analytically the exact shape of the interface of 
small particles migrating into the region initially occu- 
pied by large particles, it is assumed that the concentra- 
tion change in each section is slow enough that the fully 
radially segregated core is already completely developed. 
This approximation states a direct relationship between 
concentration C{z,t) and width and position of the re- 
gion occupied by small particles which will be given now. 

For a half-filled cylinder, the center of mass of the core 
of small particles lies below the free surface. For lower 
concentrations, we assume that all small particles occupy 
a region in the shape of a half circle, again centered at 
the same point as the full half circle. The radius r of this 
half circle is related to the concentration and the cylinder 
radius by 



The analytically obtained interfacial shape dynamics of 
the region occupied by small particles is shown in Fig. ^ 
for three different times after the start of rotation. From 
the cross sectional views of Fig. |[ we estimated that the 
front of small particles reached the left cylinder bound- 
ary around t = 45 s which gave a value of D — 0.02 
cm^/s for the diffusion coefficient in Eq. (Q). Mixing 
experiments using 6.4 mm Lucite beads gave diffusion 
coefficients around 0.1 cm^/s and the higher value is due 
to the larger particle diameter . All shapes agree well 
with the MRI data considering the level of approxima- 
tions made and the dashed line for t — 600 s predicts 
an already nearly fully radially segregated core of small 
particles which corresponds well to the right picture of 
Fig. ^. In order to derive an improved model, however, 
we need to know the exact shape of the region occupied 
by small particles and the exact position of the center of 
mass of this region. 
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FIG. 4. Theoretical prediction of the interface shape for 
three different times after start of rotation using the same 
constant diffusion coefficient D = 0.02 cm^/s. 

The proposed model with a constant diffusion coeffi- 
cient D in Eq. always leads to symmetric concentra- 
tion profiles, which can be seen from Eg. ([^). However, 
the experimental data presented in Ref. [|l5| indicate that 
even though this is a good approximation for small num- 
bers of rotations, clear deviations from symmetric profiles 
are detected for large numbers of rotations. Since parti- 
cle concentrations vary in space and time and might have 
different particle packings due to the radial segregation 
process, the mobility of the particles is also affected. The 
proposed model enables us to test the effects of a concen- 
tration dependent diffusion coefficient on the particle mo- 
tion, especially the shape dynamics. We will investigate 
two linear dependencies of the form 

D±{C) = Do ± Do{C - 0.5) (3) 

which fulfill 



(D) = r D±iC)dC^Do 
Jo 

and solve Eq. (|^) numerically by using a standard finite- 
difference procedure. 

This linear concentration dependence has a significant 
effect on the shape dynamics which is shown in Fig. |^ 
by comparing both laws with the shape given by a con- 
stant diffusion coefficient at < = 90 s. The effect is most 
pronounced in the left part of the cylinder, i.e. in the 
region the small particles propagate into. When D in- 
creases (decreases) with small particle concentration C, 
the diameter of the segregated core close to the left cylin- 
der wall is larger (smaller). Unfortunately, we are not 
able to judge, which of the investigated three concentra- 
tion dependencies of D gives the best agreement with 
the experimental data. The spatial resolution is not high 
enough to reconstruct the full shape which would be nec- 
essary in order to distinguish a dependence of the form 
D± (C) from a constant diffusion coefRcient with a larger 
(for D-) or smaller (for £)+) value. However, we will 
present in the next section a more suitable method to do 
so. 
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FIG. 5. Theoretical prediction of the interface shape at 
t = 90 s for three different concentration dependencies of the 
diffusion coefficient D. 

Since the larger particles give a higher MRI signal, 
the propagation of these particles into the region initially 
filled with small particles can be studied by recording in- 
tensity values that correspond to large particles. These 
values were extracted from the side views of Figs. || and ^ 
and are shown in Fig. |6| as filled circles with correspond- 
ing error bars. Also shown as dashed line is the theoret- 
ical curve according to Eq. (^) using the same constant 
diffusion coefRcient oi D — 0.02 cm^/s as in Fig. The 
data shows that the front advancement is underestimated 
when a constant diffusion coefficient is used. The particle 
dynamics are very different for particles in the fluidized 
surface layer and for particles following the solid body 



rotation. The fluidized surface layer has a thickness of a 
few large particles for the range of rotation speeds used 
here. For a low concentration of small particles, these can 
be trapped in a region well below this layer and they will 
show a low mobility relative to one another due to per- 
colation. For higher concentrations, the segregated core 
is large enough so that some of the small particles in the 
core are mixed with the large particles in the fluidized 
surface layer. Therefore, not all small particles can re- 
main in the segregated core but some will be recirculated 
in the flowing layer. The mobility of particles in this layer 
is much higher than for particles in the segregated core 
which makes us believe that the diffusion coefficient is 
concentration dependent. 

The exact functional dependence can only be inferred 
from a detailed understanding of the microscopic par- 
ticle motion of the two species, e.g. by looking at the 
asymmetry of the concentration profile or by calculating 
the diffusion coefficient microscopically. Unfortunately, 
neither the experimental nor the numerical results avail- 
able are accurate and detailed enough to address this 
question fully. However, it is already very instructive to 
consider the two linear dependencies proposed in Eq. (||) 
and we show in Fig. ^ as solid line the numerical result 
for D{C) = D_{C). The result for D{C) = D+{C) will 
lie between the solid and dashed lines and in either case, 
we get a better agreement with the experimental result 
when a constant diffusion coefficient D{C) = Dq is used. 
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FIG. 6. Propagation of large particles along rotation axis, 
comparison of experimental data (•) and from diffusion equa- 
tion using a constant (• • ■) or a concentration dependent ( — ) 
diffusion coefficient D. 



IV. CONCLUSIONS 

Using MRI measurements, we investigated the shape 
dynamics and interface propagation in a slowly rotating 
half-filled horizontal cylinder. By assuming that the ra- 
dial segregation and the axial migration process occur 
on different time scales, we showed that the propagating 
front can be well described by a one-dimensional diffu- 
sion process. The experimentally observed concentration 
profiles are slightly asymmetrical [ p^ , however, this can 
be understood by using a concentration dependent dif- 
fusion coefficient. We demonstrated the effect when a 
linear dependence of the diffusion coefficient on the par- 
ticle concentration is used in the analytic description. In 
this case, we found a better agreement with the experi- 
mental results as compared to the case when the diffu- 
sion coefficient is constant. This implies that the shape 
dynamics of the interface propagation can be described 
more accurately when a concentration dependent diffu- 
sion coefficient is considered. However, in order to deter- 
mine the correct concentration dependence of the diffu- 
sion coefficient, one must compare the shape of the whole 
concentration profile or one has to calculate the diffusion 
coefficient microscopically from the particle movements. 
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